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The experimental data on the ^^O+^^C and ^*0+^^C elastic scatterings and their optical model 
analysis are presented. Detailed and complete elastic angular distributions have been measured at 
the Strasbourg Vivitron accelerator at several energies covering the energy range between 5 and 10 
MeV per nucleon. The elastic scattering angular distributions show the usual diffraction pattern and 
also, at larger angles, refractive effects in the form of nuclear rainbow and associated Airy structures. 
The optical model analysis unambiguously shows the evolution of the refractive scattering pattern. 
The observed structure, namely the Airy minima, can be consistently described by a nucleus-nucleus 
potential with a deep real part and a weakly absorptive imaginary part. The difference in absorption 
in the two systems is explained by an increased imaginary (mostly surface) part of the potential 
in the ^^O+^^C system. The relation between the obtained potentials and those reported for the 
symmetrical ^^O+^^O and ^^C+^^C systems is drawn. 
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I. INTRODUCTION 



Significant progress has been made in the last decade 
in the understanding and determination of the nuclear 
optical potential for light heavy-ion systems. In the de- 
scription of heavy-ion collisions and the accompanying 
compound nucleus formation, the strongly absorptive in- 
teraction is usually present. In certain light heavy- ion 
collisions, involving closed or semiclosed shell nuclei, the 
number of open reaction channels is small and the absorp- 
tion is weaker |l| . Consequently, resonant and refractive 
phenomena have been observed in these systems. Recent 
reviews of theoretical and experimental results on res- 
onant and refractive phenomena can be found in 
Reported systematic measurements of the elastic scat- 
tering which cover broad angular and energy ranges and 
their optical model analysis have resulted in a deeper 
insight into the dynamics of light heavy-ion scattering. 
Such progress is closely connected with the observation 
of refractive effects. It has been shown that the presence 
of the nuclear rainbow and, the associated Airy struc- 
ture in elastic scattering reduces the ambiguities of the 
optical potential character (shallow or deep). The exis- 
tence of systematic studies, where the change of potential 
parameters can be followed as a function of energy, is es- 
sential for the determination of a unique nucleus- nucleus 
potential at different energies and in different systems. 



The main features of the optical rainbow lie in the re- 
fraction and reflection. In the case of interference due to 
the contributions around the maximum deflection angle, 
higher-order maxima, i.e.. Airy structure, appear inside 
the lighted region. The observation of the rainbow pat- 
tern in heavy-ion collisions is closely connected with the 
degree of transparency. Water is transparent to visible 
light, while the nuclear rainbow is damped by the pres- 
ence of absorption. The rainbow scattering appears when 
the nuclear potential is strong enough to deflect particles 
into negative angles (so-called farside) . The combination 
of these two features, deep real potential and incomplete 
absorption, makes possible the observation of distinct re- 
fractive effects like the Airy minima, the rainbow angle 
and, the rainbow dark side. 

The energy range between 5 and 10 MeV per nucleon 
can be defined as an intermediate region lying between 
that where the molecular-resonant states show up and 
that with the prominent appearance of the nuclear rain- 
bow. The experimental results on binary channels re- 
ported up to very recently in the literature for this en- 
ergy range are rather scarce. The refractive studies have 
mainly involved symmetrical systems, mostly ^^C-l-^'^C 
and 160-^160 To cover the lack of system- 

atic measurements, several experiments have been per- 
formed at the Strasbourg Vivitron accelerator, which is 
well suited for studies in this intermediate energy range. 
Detailed and complete elastic angular distributions have 



1 



been measured for three systems: "'^^0+"'^^0, ""^^O+^'^^C, 
and In this work we concentrate on the study 

of the nonidentical systems where the angles beyond 90° 
can be explored directly. 

In the literature, the data for the system are 

not as complete and numerous as for identical boson sys- 
tems. Relative to the ^^O beam energies of the present 
work, the elastic scattering angular distributions have 
been reported and analyzed at higher energies: i?iab — 
132 ll, 139, 216, and 311 MeV |lj, 608 @, and 1503 
MeV ]lj, and at lower energies: i?iab = 24, 35, 45 
65, and 80 ||l|,|l3 MeV. The data reported in Refs. 
5||l6| cover only limited angular ranges and do not 
reveal any significant refractive features. The data at 608 
and 1503 MeV have been shown to contain some refrac- 
tive effects displaying sensitivity to the features of the 
real potential. Recently, the elastic scattering measure- 
ments over a wide angular range have been published at 
energies just above the energy range covered in our study, 
i.e. at i?iab=132 (re-measured and extended to large an- 
gles), 170, 200, and 230 MeV |l|. No refractive study of 
the ^^0-|-^^C system has been reported in the literature 
so far. 

Our elastic scattering data and the optical model anal- 
ysis for the ^^0-1-^^0 symmetric system were presented 
in |l|-0, and for the system in |2|. The pre- 

liminary results on the ^^0-|-^^C scattering we presented 
in In this article, we report on an additional 

optical potential analysis of the ^'^O-I-^^C and ^^O-I-^^C 
systems. In particular, we present the nearside- farside 
decomposition of the deduced potentials which allowed 
one to identify Airy minima and to follow their evolution 
as a function of beam energy. To have a better connec- 
tion with the higher-energy domain, in the analysis of the 
^^0-|-^^C system we have included the lowest energy (132 
MeV) elastic scattering angular distribution reported in 
. The change of the potential due to the addition of 
two extra neutrons in the ^^0-|-^^C system will be explic- 
itly discussed. 



—70°) at 7.8 cm from the target, a self-supporting 
carbon film - 20//g/cm2 thick. The and 
^*0-|-^^C elastic scatterings were measured at beam en- 
ergies £;iab(^^0)=62, 75, 80, 94.8, 100, 115.9, and 124 
MeV and, at £:iab(^*0)=66.2, 85, 100, and 120 MeV, re- 
spectively. The electronics and data acquisition system 
of the Q3D spectrometer with the proportional counter in 
its focal plane and of the PS detectors in the kinematical- 
coincidence mode were independent of each other. The 
two detection systems had an overlapping angular re- 
gion. The obtained center-of-mass elastic scattering an- 
gular distributions span the angular range between 10° 
and more than 140° for both systems studied with good 
statistics at all measured energies. 

All angular distributions and the optical model de- 
scription of the system were reported in [ p2| . 
In the present work we intend to study in more detail the 
connection with the higher-energy region. Thus, we are 
not going to discuss the results obtained at the two low- 
est energies, 62 and 75 MeV. In fact, it has recently been 
shown that at even lower energies the deep real potential 
describes both the observed structure in the fusion exci- 
tation function and the gross features of elastic angular 
distributions |p7|| . 

Resonant phenomena in heavy-ion reactions are usu- 
ally restricted to systems composed of the so-called a- 
particle nuclei. Considering the presence of resonances 
as a signature of weak absorption, one would expect the 
refractive effects to appear in the higher-energy domain 
for the same systems. A description of all details that 
appear in angular distributions where the resonant ef- 
fects are superimposed on the gross structure is unlikely 
to be achieved with a mean field approach, i.e., the opti- 
cal model. This is the case with the ^^0-|-^^C system at 
£^iab=66.12 MeV, the lowest energy of the present mea- 
surement. At this energy a broad resonance has been re- 
ported and the energy spectra and angular distributions 
of the binary channels were discussed in [Ea,E7| . 



II. EXPERIMENTS 

The three reactions ^^O+'^^C, and ^^O+^^O 

were studied under the same experimental condition, 
which will be briefly descri bed in this section. More de- 
tails can be found in Refs. 
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The elastic scattering angular distributions were mea- 
sured at laboratory angles between 5° and 20° in steps 
of 0.5 degree using the Q3D magnetic spectrometer 
and its associated detection system. At larger angles 
all binary channels were recorded simultaneously using 
a fixed kinematical-coincidence setup composed of two 
position-sensitive (PS) silicon detectors. The setup al- 
lowed complete mass and Q-value identification of the 
binary fragments. PS detectors were placed on both 
sides of the beam (15° < 6* < 50° and -35° < 9 < 



III. OPTICAL MODEL POTENTIAL 

In the analysis of the ^^O+^^O and 
p2t systems several choices were tried for the real part 
of the potential. Our understanding of the proper- 
ties and microscopic grounds of the optical potential is 
largely based upon the folding model using a realistic 
nucleon-nucleon effective interaction The energy- 

and density-dependent interaction was used to gener- 
ate a microscopic potential, which together with a phe- 
nomenological imaginary term, successfully describes our 
i6o_|_i6o and ^'^O-I-^^C data as well as the collisions in 
slightly higher and lower energy domains 1^,0. The 
resulting folded potentials are strongly attractive, and 
weakly energy dependent. It has been shown (see Fig. 1 
in [pll) that the square of the Woods-Saxon form factor 
(WS2) with a suitable choice of radius and diffusencss 
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gives a shape that is very close to the folded potential. 
Indeed, equivalent fits to the data are found with either 
the folded or the WS2 potentials. In the present case, 
the real part of the potential is taken in the WS2 form: 



V{r) - -V[fv{r)f 



/y = [1 + exp{- 



(1) 



Rv . 



av 

The shape of the folded potential for the system 
is closely reproduced by Eq. (|l|) with the radius i?v =4.0 
fm and the surface diffuseness ay = 1-4 fm. These values 
were kept constant throughout the best-fit automatic pa- 
rameter searching procedure. In the energy range consid- 
ered, no significant improvement of the fit was obtained 
by allowing i?v and ay to vary freely. Slightly different 
values of the radius and diffuseness of the real potential 
for ^^0-|-^^C were adopted owing to the addition of the 
two neutrons. They were kept constant throughout the 
fitting procedures at the values Ry =4.08 fm and ay = 
1.38 fm. 

Two choices were used for the imaginary potential, 
a pure "volume" term of the Woods-Saxon (WS) type 
(VFws) and a volume term together with an additional 
"surface" term (Wws2 + W^d)- For the second choice of 
the imaginary potential we adopted the WS2 form fac- 
tor for the volume term, and the radial derivate of the 
WS form for the additional surface term. The imaginary 
potentials are thus given by the following equations: 



W{r) = W^s{r) = -iW[Jw{r)l 
W{r)=W^f,s2{r) + WT,{r) 



(2) 



(3) 
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In the following text, the optical potentials defined by 
Eqs. (0) and (||) will be called parameterization P(Wws) 
and those defined by Eqs. (|l|) and (H) will be called 
parameterization P{Wws2 + W^d)- The fits were ob- 
tained through the automatic search option in the pro- 
gram PTOLEMY ||. 

The usual criterion was used to judge the quality 
of agreement with the data: 
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(4) 



The (Tth and <7ex are theoretical and experimental cross 
sections, Auex are the uncertainties in the experimental 
cross sections, is the total number of angles at which 



measurements were carried out, and iVp is the number 
of free fitting parameters. Two choices for the cross sec- 
tion uncertainties were considered. We either used the 
experimental cross section uncertainties or a fixed per- 
centage of the experimental cross section for all measured 
angles. Compared with forward angles, the experimen- 
tal cross-section uncertainties can be relatively large at 
intermediate and backward angles owing to the smaller 
cross section. Thus, the searching procedure does not 
give sufficient weight to those angles at which the refrac- 
tive effects are expected to be observed. We have found 
that using a uniform percentage (10%) a better fit at 
intermediate and backward angles can be obtained. 

For an interaction potential lJ-^(r) between nuclei with 
nucleon numbers A\ and A2, the volume integral 



in 
A1A2 



U-E{r)r'^ dr 



(5) 



is a sensitive measure of the potential strength. This def- 
inition applies to the real and imaginary parts of Usir). 
Generally, the optical model analyses can give several 
discrete families of real potentials that fit an angular dis- 
tribution equally well. It has been shown that the volume 
integrals are valuable tools for classifying different optical 
model families, especially in terms of global systematics, 
when different systems are compared. 

The simple optical model potential is a local, angular- 
momentum independent potential. Thus, effects like the 
elastic transfer, an a-transfer in or a ^He- 

transfer in ^^0-|-^^C, which can introduce additional 
structure in the angular distribution, are not included. 

The obtained angular distributions are complicated to 
interpret owing to the contributions and interferences 
of many involved partial waves. For this purpose, the 
nearside-farside decomposition technique, first presented 
by Fuller | |30| ], is a helpful method for interpreting the 
optical model results. The nearside-farside decomposi- 
tion separates the trajectories which originate from two 
different sides of the scattering potential. By convention, 
the nearside and farside correspond to classical trajec- 
tories with positive and negative deflection angles, re- 
spectively. The interference between nearside and far- 
side trajectories will lead to the Fraunhofer diffraction 
pattern, and the interference between two farside com- 
ponents (the one more peripheral and the other more in 
the interior) will lead to the rainbow pattern. The refrac- 
tive effects, namely the Airy minima and maxima, will 
therefore be present in the farside component. The cru- 
cial condition for the observation of refractive effects is 
incomplete absorption. The absorption should be weak 
enough ("transparency") to allow the inner farside com- 
ponent to be effective. The appearance of Airy oscilla- 
tions in an angular distribution is a direct indication of 
interferences between waves which are bent by the nu- 
clear mean field and provides a unique information on 
the potential at small interaction radii. 
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IV. ELASTIC ANGULAR DISTRIBUTIONS 
A. Optical model description 

Figures |l| and || show the data and the phenomenolog- 
ical model fits (thick-solid curves) of the angular distri- 
butions for the and systems, respec- 
tively. The parameters obtained using the fitting proce- 
dure are listed in Table | (parameters P(Wws2 + Wd) 
which are used on the left panels of both figures) and 
Table || (P(WVs) and the right panels of Figs. and|). 

The main features of the measured angular distribu- 
tions, namely, the forward fine Fraunhofer diffractive os- 
cillations and the broad structure at larger angles are 
well described by the calculation. An increase in cross 
sections at large angles is predicted by the calculations, 
but is smaller in some cases than observed experimen- 
tally. 

In our previous analyses of the ^^0-1-^^0 system Q, it 
was shown that for the data at the three highest energies 
(124, 115.9, and 103.1 MeV) an acceptable description 
was possible using just a volume term for the imaginary 
part. The fits at lower energies (lower than 100 MeV) 
required inclusion of a surface imaginary term. To be 
consistent, in the present analyses, the ^^O+^^O data at 
all energies were analyzed using both imaginary terms, 
volume and surface. The resulting imaginary potential is 
mostly of the volume type at higher energies. At lower 
energies the imaginary potential has a sharp-edged vol- 
ume term and a surface term which peaks outside the 
volume one (see the top-right panel of Fig. 0) . The ini- 
tial set of parameters for the fits of the system 
was the one which gives the same volume integrals of the 
real and imaginary parts (Eq. ^) as those obtained in 
the fitting procedure for the ^^0-1-^^0 system. The final 
difference between the volume integrals of the -'^^0-1-^^0 
and systems at 124 MeV is -10%. In the op- 

tical model analysis of the ^^0-|-^^C system presented in 
p2t , at all measured energies the folding and WS2 terms 
were used for the real potential and the sum of the vol- 
ume (WS) and surface terms for the imaginary potential. 
It has to be noted that the imaginary volume term in |2|] 
was of the WS type, while in the present analysis it is of 
the WS2 type when the surface term is included, which 
results in slightly different parameters. The main goal 
of this work is to establish a connection with the higher- 
energy domain and to identify the Airy minima and their 
order using the nearside-farside decomposition technique. 
To fulfill this task, we included in our study the angular 
distribution at 132 MeV ||l^. The optical model analysis 
presented in was performed by using for the real part 
either the folding or the phenomenological Woods-Saxon 
potentials. For the imaginary part, only a volume term 
was used (WS at all energies and WS2 at several ener- 
gies). Such a potential, with only a volume imaginary 
term, describes the main features of the measured distri- 
butions, but fails to describe the backward- angle oscilla- 



tions which become more important as energy decreases. 
Such a potential underestimates the observed increase of 
the yield at backward angles. 

Actually, a satisfactory description of the whole struc- 
ture observed in the ^^0-|-^^C system in the energy range 
considered, required the inclusion of an imaginary sur- 
face term. To justify this statement, the analysis with 
and without the imaginary surface term is consistently 
carried out and presented. Indeed, the distributions cal- 
culated without the surface term resulted in too little 
structure at large angles (Figs. ^ and ||, right panels). 
In addition, the predicted increase of yield at backward 
angles (in agreement with data) is less pronounced when 
the cross section is calculated using the imaginary vol- 
ume term only. It seems that for a certain energy range, 
the data cannot be sufficiently well described without 
the surface term and that the ^''O+^^C and ^^O-I-^^C 
systems at the energies considered here are within this 
range. 

For the ^^O-H^^C system, the fits were obtained start- 
ing from the parameters of the ^^0-|-^^C system. The 
results of this procedure at an energy of 100 MeV, can 
be seen in Fig. ||. The top panel shows the best fit of 
the -'^^0-1- "'^^C system with two different sets of parame- 
ters, using the sum of the volume and surface terms (solid 
curve) and using the volume term only (dashed curve). 
These sets of optical potential parameters were used to 
calculate the angular dependence for the ^^C sys- 

tem. The results of the calculation are superiniposed on 
the data in the middle panel of Fig. g. Agree- 

ment between measured and calculated values is reason- 
able at forward angles, but the calculated angular dis- 
tributions cannot describe the two orders of magnitude 
smaller yield in the ^*0-|-"'^^C system at larger angles. It 
is worth noting that even without fitting, the number of 
oscillations in the calculation obtained using the imag- 
inary term composed of the volume and surface terms 
(solid curve) agrees with the measurement throughout 
the angular range and over eight orders of magnitude of 
the differential cross section. The bottom panel of Fig. 
^ shows the results of the fitting procedure. The dif- 
ferences between the potentials of the two systems are 
presented in Fig. ^: The left-bottom panel refers to the 
P(Wws) parametrization and the right-bottom panel to 
the P{Wws2 + Wu) one. The real part remains almost 
unchanged, but the imaginary volume integral increases 
substantially from ^^^O+^^C to To describe the 

decrease of the yield at backward angles, the absorption 
has to be increased. In the imaginary part of the poten- 
tial which is composed of the volume and surface terms, 
the surface term now becomes dominant at large inter- 
action radii. 

The obtained potentials have deep real parts. When 
the additional surface term in the imaginary potential 
is included, the volume term has a smaller radius than 
the surface term and as a consequence, the volume imag- 
inary potential has an almost square profile whose de- 
tailed shape depends on the system and energy. 
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The volume diffuseness of the imaginary potential 
tends to be very small when the surface imaginary term is 
included in the potential. The volume imaginary diffuse- 
ness could be fixed at some small values, such as 0.1 fm, 
without any significant change of the fits, but we rather 
report the values resulting from the automatic search. 

The importance of the volume integrals can be shown 
by comparing the results of the two different parame- 
terizations, P(Wws) and P(Wws2 + Wb)- The poten- 
tials which belong to P(VFws) and P(Wws2 + W^) at 
100 MeV are shown in Fig. ^: left-top panel for the 
i6o_|_i2c and right-top panel for the ^^O-I-^^C systems, 
respectively. The dashed curves correspond to P{Wws): 
while the solid curves correspond to the P(Wws2 + W^d)- 
The two solutions give the same (within a few percent) 
volume integrals (see Table pT[ ) . In the imaginary part, 
when the individual surface and volume terms are con- 
sidered, both terms are important. As mentioned above, 
the common characteristic of the Wws2 + Wu solutions 
is that the volume imaginary term tends to have a sharp 
edge and a smaller radius than the surface term. The pre- 
cise balance between the two imaginary terms is not so 
well established as the values of the integrals. The imag- 
inary potential consisting of the volume term only, tends 
to smooth out the potential pocket of the P(Wws2 + Wd) 
in a way as to keep the imaginary volume integral about 
the same. 

The volume integrals of the real potentials are very 
similar for both parameterizations and also for the phe- 
nomenological and the folding potentials. In |jl^ two 
families of discrete WS sets were obtained using the fit- 
ting procedure. The obtained real WS potentials have 
different depths and similar radii and consequently dif- 
ferent volume integrals. At lower energies both WS sets 
describe the data equally well, but the deeper potentials 
(the WS2 set in iQ) fail to reproduce the right order of 
Airy minima at higher energies. This is clearly visible in 
the data at 200 MeV of |1^ where the first Airy mini- 
mum and the rainbow angle are observed. Comparing the 
values of the volume integrals and the numerical values 
of the potentials, we may conclude that the real poten- 
tial obtained in this work belongs to the same family (the 
WSl set in [Q). The selection of a unique potential fam- 
ily results in identification of the order of Airy minima. 
Furthermore, it makes possible to follow the evolution of 
the Airy structure with energy. 



B. Nearside-farside decomposition 

The nearside-farside decomposition is particularly 
helpful in the intermediate energy range of our study. 
In Figs. and ^ the thin-solid and dashed curves show 
the nearside-farside decompositions of the obtained opti- 
cal model fits of the ^'^O-f ^^C and the ^^O-I-^^C systems, 
respectively. 

The angular distributions at the highest and the lowest 



energy of the ^^0-|-^^C system (top and bottom panels 
of Fig. 1^) show different patterns. The fit and calcu- 
lation for the 132 MeV angular distribution display the 
features which are common to this and similar systems at 
higher energy. We can describe the angular distribution 
by Fraunhofer oscillations at forward angles, glory effects 
at backward angles, and intermediate structures, namely 
the deep minimum at 80°, as an Airy oscillation. The 
nearside-farside crossover takes place around 25°. Be- 
yond this angle, the angular distribution is farside domi- 
nated and the structure of the distribution is refractive in 
its origin, i.e., it is the result of the interference between 
the two farside subamplitudes. By definition, the near- 
side and farside amplitudes are equal at 180° and their 
interference produces strong and rapid oscillations (glory 
effect) in the angular distributions as we approach the 
largest angles. The situation changes at lower energies. 
The nearside-farside crossover moves to larger angles and 
the Fraunhofer oscillations also extend up to large angles, 
they are superimposed on broad oscillations of refractive 
origin. 

The comparison between the decompositions of the two 
different sets of parameters of the imaginary potential 
(left and right panels of Fig. |l|) shows that the obtained 
farside components are very similar, while the nearside 
components exhibit quite different behavior. As already 
stressed, the refractive effects appear in the farside com- 
ponents and, omitting the different oscillatory pattern at 
very backward angles, the location of the Airy minima 
of the two different parameterizations remains the same. 
Such behavior clearly shows that the Airy structure is 
determined entirely by the real potential. 

The nearside and farside components are of compara- 
ble strengths at intermediate angles for the P{WwS2 + 
Wu)- The nearside recovers its strength beyond the 
nearside-farside crossover. For the decomposition of the 
calculations obtained with the P{Wws) and shown on 
the right panels, the farside components behave more as 
was observed at higher energies, i.e., the angular distri- 
bution is farside dominated, whereas the nearside is more 
than two orders of magnitude smaller. 

The angular distribution measured at 115.9 MeV is a 
typical example of the cross section in this intermedi- 
ate energy range. The nearside-farside crossover takes 
place at 50°, 60°, and 100° for the P{Wws2 + Wu)- 
The nearside and farside amplitudes become interlaced 
owing to their similar strengths and intersect at several 
points. The nearside-farside crossover takes place at 27° 
for the P(Wws) and above this angle the distribution 
is farside dominated. Again at backward angles, the in- 
terference between the nearside (coming partly from the 
farside component after passing through 180° and going 
around) and farside components is clearly visible in the 
angular distribution. Looking only at the farside ampli- 
tude, two Airy minima appear in the amplitude {A2 near 
95° and A3 near 60°) for both parameterizations. The 
data themselves show a wide minimumlike structure at 
these angles. 
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Owing to the stronger absorption in the 
system (see Fig. H), even at our highest energy (120 
MeV) , the nearside and farside amphtudes are of similar 
strength throughout the whole measured angular range 
(top- left panel in Fig. ||) for the P{Wy/s2 + Wu), whereas 
for the P{Wws) the nearside-farside crossover is shifted 
toward larger angles (top-right panel in Fig. ^). The 
similarity of the nearside and farside strengths prolon- 
gates the appearance of the diffractivelike oscillation up 
to large angles. Moreover, the nearside-farside crossover 
moves even to larger angles as energy decreases. In the 
description of these data we have to keep in mind that 
the farside which turns over at 180° becomes the near- 
side, and that the structure observed at backward angles 
can be refractive in its origin [^l|j22| ]. In spite of this 
angle-extended Fraunhofer diffraction pattern the Airy 
minima can be discerned by considering only the farside 
amplitudes. In the raw data themselves the refractive 
effects are observed rather as irregularities in the oscil- 
lation pattern. Owing to the stronger absorption in the 
^*0-|-^^C system, the nucleus-nucleus interaction trans- 
parency is reduced and the interference between the two 
farside subamplitudes does not dominate the angular dis- 
tributions in the energy range studied. The nearside- 
farside decomposition with the P(WwS2 + Wu) shows 
that the experimental distributions are nearside domi- 
nated, i.e., the diffraction is important throughout the 
measured angular range. We have to keep in mind that 
the P(H^ws2 + Wu) better describes the observed struc- 
ture at larger angles. Nevertheless, even if the deep min- 
ima of refractive origin in the angular distribution are 
not observed, the description of the structure is within 
the refractive picture, i.e., the deep real potential is re- 
quired. We want to emphasize that whenever a choice 
among alternative solutions of the optical potential had 
to be made, the main criterion was the regularity (mono- 
tonic variation) of the parameters. The regularity was 
based on the expectation that the potential should not 
change rapidly with energy and for neighboring systems. 
The resulting potential parameters which show such a 
regularity allow the connection of the observed structure 
in the elastic angular distributions of the ^^0-|-^^C sys- 
tem with the refractive effects, as observed in the neigh- 
boring and more "transparent" ^^0-1-^^0 and ^^0-|-^^C 
systems. Since the refractive effects are not so prominent 
in the ^*0-|-^^C angular distributions, the optical model 
ambiguities remain to a certain extent. 

We may conclude that in the intermediate energy 
range, between 5 and 10 MeV per nucleon, refractive ef- 
fects are an important ingredient of the global structure 
of an angular distribution. The ambiguities in the deter- 
mination of these refractive effects are closely related to 
the degree of absorption in the system. Even if the re- 
fractive effects are masked, the structure observed in the 
measured angular distributions can be at least partially 
explained as being of refractive origin. The stronger the 
absorption, the more masked are the refractive effects 
and their description becomes more complicated and less 



accurate. Alternative explanations of the observed struc- 
tures in angular distributions at this intermediate energy 
range can be helpful. 

Another possibility is the decomposition into the 
barrier- and internal- wave (B/I) components [Bl 32 1. The 



B /I decomposition makes sense if the real part of the po- 
tential is deep enough for the effective potentials to dis- 
play a " potential pocket" and if the absorption is incom- 
plete. Recently, the B/I decomposition has been applied 
to our data, at all energies for ^^O+^^O |^,^, and also 
at one representative energy for our nonidcntical systems 
16q_^12c j^j^^ iSQ-^i^c Q. Within the B/I decomposi- 
tion, the Airy minima are explained as an interference of 
the barrier-wave and internal-wave subamplitudes. The 
results of the B/I decompositions for the lower and higher 
energy range of our ^^0-1-^^0 data differ only in the rela- 
tive importance of the two components. Inclusion of the 
surface imaginary term modifies the internal contribu- 
tion only slightly, but produces the increase of the barrier 
cross section at large angles, which results in the increase 
of the yield and the presence of oscillations in the cross 
section. The prominence of the Airy minimum in the 132 
MeV angular distribution of ^^0-|-^^C is due to similar 
magnitudes of the internal-wave and barrier-wave com- 
ponents in the angular range around 80°. Owing to the 
stronger absorption in the ^*0-|-"'^^C system, the internal- 
wave cross section is about three orders of magnitude 
lower than in the ^^0+^^C system, which almost com- 
pletely smears out the appearance of Airy minima in the 
experimental angular distributions. 

Very recently a semiclassical analysis of the Airy-like 
pattern has been presented [Q. The oscillations in the 
farside subamplitude are explained as a result of inter- 
ference of the first and second term of a multireflection 
expansion of the scattering function. The physical con- 
tents and the obtained decomposition are very similar to 
the results of the B/I decomposition. 



C. Interpretation of the results in terms of Airy 
minima 



The rainbow and associated Airy oscillations in the 
scattering appear in the farside component. The Airy 
structure is controlled by the real potential and its ap- 
pearance is often obscured by the presence of absorption. 
As the energy increases, the rainbow and Airy minima 
move forward in angle. The identification of the struc- 
tures of refractive origin can be simplified by using a re- 
duced imaginary potential. Figures ^ and |6| show the far- 
side amplitudes with a 50% reduced imaginary strength 
of the ^^O+'^'^C and ^^O-I-^^C systems, respectively. The 
real potential is weakly dependent on energy and the 
smooth energy evolution of the Airy minima can be ob- 
served in the calculated distributions in Figs. || and ||. 

Let us first discuss the results for the ^^0-|-^^C system 
(Fig. |). In Refs. ||l|,||| the strong minimum at 80° in 
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the measured 132 MeV angular distribution is identified 
as the second Airy minimum A2. The minimum at 55° 
was identified as A3 and that at 35° as A4. As already 
discussed, the real potentials obtained in our analysis for 
the system belong to the same family as the real 

potentials reported in . Following the systematics of 
the order of Airy minima, the minima in the angular dis- 
tribution at 124 MeV appearing at angles 88°, 60°, and 
40° are identified as A2, A3 and A4. At 100 MeV the 
A2 minimum moves to 120°, A3 to 80°, and A4 to 50°. 
At very forward angles, around 30°, the fifth Airy mini- 
mum can also be discerned in the farside component at 
100 MeV. Of course, it is not easy to see the remnants 
of all these minima in the actual data, especially for the 
higher-order minima appearing at more forward angles 
where strong Fraunhofer oscillations are dominant. 

No detailed elastic angular distributions at the 
^*0-|-^^C system for energies higher than those in the 
present study (where one expects a more favorable situ- 
ation for the identification of the Airy minima and their 
order) have been reported so far. In the identification of 
the order of Airy minima, we assume a similar behavior 
in the and systems. Therefore, it is 

likely that the A2 and A3 minima appear in 
at energies around 100 MeV as well. According to this 
assumption, the minimum around 105° in the 120 MeV 
angular distribution is labeled as A2, and the minima at 
75°, 50°, and at the very forward angle of 30° as A3, A4, 
and (see Fig. The minima move to larger angles 
at lower energies and at 100 MeV the A3 moves to 105°, 
A4 to 72°, and A5 to 45°. At 85 MeV, even the Aq can 
be discerned at a very forward angle. 

The positions of the Airy minima as a function of 
center-of-mass energy for both systems studied are shown 
m Figs. and |. A very regular parabolic behavior is 
observed. 

The angles at which the minima take place in the far- 
side amplitudes do not depend on the imaginary poten- 
tial strength and no important angular shift is observed if 
the absorption is switched on. The depth of the minima 
in the farside changes with absorption. The absorption 
dictates which minimum will be dominant in the farside 
amplitude. The stronger the absorption, the deeper the 
minima will be at larger angles and conversely at smaller 
angles. It has to be noted that the interference between 
the volume and surface components introduces some ad- 
ditional small oscillations around the Airy minima. This 
speaks in favor of the conclusion that the volume and 
surface absorptive components play different roles. 



D. Potentials, volume integrals, and global 
systematics 



In this section we present all optical potentials ob- 
tained for the four systems considered in our study: 
I'^O+ieo ll9|-|2ll], 160+12C [p2|-p4i, i^O+i^C MM, 



and 12C-I-12C [@,|4|j3|. Figure | shows the energy and 
system variations of the obtained potentials (real and 
imaginary parts) at representative energies of the energy 
range studied. 

In all systems, a deep real part is required for a good 
description of the elastic angular distributions. The ob- 
served structure in the distributions is explained through 
refractive effects. The imaginary part changes from sys- 
tem to system. While the inclusion of the surface term 
was needed only for the lower energy range in ^^O+^^O, 
the ^^0-|-^^C system required a surface term for all en- 
ergies, and without it the structure at larger angles can- 
not be well described. The surface term becomes even 
more important in the obtained imaginary potential of 
the ^*0-|-^^C system. When the surface term is needed, 
the absorptive potential has a characteristic shape: a vol- 
ume term which tends to have a small diffuseness plus a 
Gaussian-like surface peak. 

The most typical examples of such a shape are the op- 
tical potentials for ^^O+^^O at 75 MeV, for 
at 124 MeV, and for ^^O-I-^^C at all energies studied. 
Within a simple interpretation, we can associate the vol- 
ume term with the absorption due to fusion and related 
processes, and the surface term with more direct reac- 
tions. These different shapes of the absorptive potential 
can be explained as a consequence of the evolution of the 
reaction mechanisms which change from the fusion type 
to a more direct type of reaction. This enhancement of 
the surface integral is a signature of a larger number of 
open direct reaction channels and is in agreement with 
the number of open channel calculations Comple- 
mentary studies of the distribution of the incident flux 
into the different available exit binary channels, as well 
as of the underlying reaction mechanisms are actually 
underway for the ^^0-|-^^C and ^^0-|-^^C reactions. 

Tabic III lists the real and imaginary volume integrals 
of the potentials obtained for the oxygen-f carbon sys- 
tems studied and for the two adopted sets of parame- 
ters. The real and imaginary volume integrals display a 
smooth behavior with energy. The real volume integral 
generally decreases with energy, while the imaginary vol- 
ume integral increases. Such behavior is in agreement 
with the dispersion relation predictions. Although the 
shapes of the imaginary potential are radically different 
for the two sets of parameters, the full imaginary volume 



integrals (columns 5 and 7 of Table [II) are about the 
same. 

In the review article of Brandan and Satchlcr Q , the 
values of the real and imaginary volume integrals as a 
function of projectile energy per nucleon for the poten- 
tials that fit the i^C-^i^C, 160-^160, and data 
have been extracted (see Fig. 6.7 of [Q). It is shown 
that the different real potentials all have similar volume 
integrals and that the imaginary volume integrals show 
a smooth behavior as a function of energy in spite of 
various parameterizations adopted for the imaginary po- 
tentials. It is quite remarkable to notice that the volume 
integrals reported in our work closely follow the global 
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systematics evidenced in 



V. CONCLUSION 

In this paper we have presented the elastic scatter- 
ing data of the ^^O+^^c and ^^O+^^C reactions, to- 
gether with the optical model analysis. We have used 
the nearside-farside decomposition technique to inter- 
pret the complicated features of the angular distribu- 
tions. The elastic angular distributions show not only the 
usual Fraunhofer diffraction pattern, but also, at larger 
angles, refractive effects in the form of nuclear rainbow 
Airy structures. The main features of the measurements 
are very well described by the optical model fits. A deep 
real part is required for the description of refractive ef- 
fects. The imaginary potential is weakly absorptive and 
reflects the presence of incomplete absorption. The two 
studied systems, ^^O-f^^C and ^^0-|-^^C, have different 
shapes of the imaginary potentials and different values of 
the volume imaginary integrals. The volume imaginary 
integral increases for the ^*0-|-^^C system, mostly in the 
imaginary surface part. The enhancement of the surface 
integral is a signature of a larger number of open direct 
reaction channels. This conclusion is supported by the 
study of the inelastic and transfer channels, as well as by 
the number of open channel calculations. 

The inclusion of the surface imaginary term provides 
the needed increase of the yield at large angles, but also 
defines a characteristic shape of the imaginary potential. 
When the surface term is required by the data, the imagi- 
nary volume term tends to have a small difusseness and a 
smaller radius than the surface imaginary term. One can 
imagine that the reflection from a potential with such a 
profile may produce additional interferences resulting in 
more structured angular distributions at larger angles. 
As has been mentioned, the effects of elastic transfer 
could also cause an increase of the yield at the largest an- 
gles and additional structure in the angular distributions. 
The specific imaginary potential obtained here can be un- 
derstood through coupling effects which should have an 
impact on the elastic angular distribution. An additional 
analysis like the explicit inclusion of elastic transfer could 
provide a better understanding of the origin of the large 
angle structures, but it is unlikely that this would appre- 
ciably change the results for the real part of the poten- 
tial. Therefore, the most important conclusion about the 
depth of the real potential, the potential which describes 
the refractive effects and defines the position of the Airy 
minima, will remain unchanged. 

The appearance of refractive structures permits an 
unambiguous determination of the optical potential. 
Within this analysis we do not claim that the obtained 
potentials are unique, energy by energy. What we believe 
as unique is the regularity and systematics of the poten- 
tials as a function of energy and also of the target and 
projectile matter distribution. The obtained potentials 



for the light heavy-ion systems studied have a deep real 
part and a weakly absorbing imaginary part. In spite of 
the differences in the shape of the imaginary potentials 
for the i^C-^i^C, ^^O+^^O, ^^O-f ^^C, and sys- 
tems (see Fig. ^ the volume imaginary integrals agree 
with the global systematics (evidenced in Q). The sys- 
tematics is based on the refractive nature of light heavy- 
ion collisions and is in agreement with the dispersion re- 
lation predictions. The obtained volume integrals, both 
real and imaginary, agree with potentials which fit the 
higher and lower energy data. 

We may thus conclude that in the intermediate energy 
region, between 5 and 10 MeV per nucleon, the refrac- 
tive effects, although masked partially by diffraction, are 
still strong enough to permit a simple mean-field optical- 
model explanation of the main structures of the elastic 
scattering cross sections. 

VI. ACKNOWLEDGMENTS 

We are indebted to A. A. Ogloblin and his group for 
making their experimental angular distributions avail- 
able to us. We would also like to thank Dao T. Khoa 
for providing us with his nearside/farside decomposition 
codes as well as with a UNIX adapted version of the code 
PTOLEMY. 



8 



TABLE I. Phenomenological potentials, the real part is a 
WS2 term and the imaginary part is a sum of the WS2 (vol- 
ume) plus WSD (surface) term. V and the subscript V stand 
for the real term, W and the subscript W for the volume 
imaginary term and the subscript D for the surface imagi- 
nary term. The aw could be fixed arbitrarily at some small 
value, such as 0.1 fm, without much influence on the fits, but 
we prefer to report the values resulting from the automatic 
search. 
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TABLE IIL Volume integrals for the potentials obtained in 
i6q_^i2(-, 180-^12^. Wws2 + Wd stands for the param- 
eters listed in Table ^ and Wws for the parameters in Table 
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74 
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TABLE II. Phenomenological potentials, the real part is 
a WS2 term and the imaginary part is a WSl term (pure 
volume). 
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FIG. 1. Elastic scattering data shown as the ratio to 
the Rutherford cross sections and optical model calculations 
(thick-solid curves) with phenomenological potentials (left 
panels with parameters from Table |, and right panels from 
Table 0) of the ^^O+^^C at 132, 124, 115.9, and 100 MeV 
(from top to bottom, left and right, respectively). The near- 
side (thin-solid curve) and farside (dashed curve) subampli- 
tudes for different potentials are plotted. 

FIG. 2. Same caption as for Fig. |l| but for the ^*0-t-^^C 
elastic scattering at 120, 100, and 85 MeV. 



FIG. 4. The real and imaginary potentials at 100 MeV of 
the ^''O-f ^^C (left-top panel, and solid curves in bottom pan- 
els) and ^*0-|-^^C (right-top panel and dashed curves in bot- 
tom panels) for the two sets of parameters listed in Tables | 
andy. 

FIG. 5. The calculated farside amplitudes for all measured 
energies of the ^^0-|-^^C system using the parameters from 
Table |. The imaginary strength has been reduced by 50% to 
emphasize the refractive effects. Curves at different energies 
have been shifted by a factor of 100 for clarity. 



FIG. 3. Measured elastic angular distributions displayed 
as ratio to the Rutherford scattering of the ^'^O-I-'^^C (stars, 
panel a) and '^^O+^^C scattering (dots, panels b and c) at 
an incident energy of 100 MeV. The solid curve (parameters 
from Table ^ and dashed curve (parameters from Table ^ 
of panels a and c represent optical-potential fits. The same 
fits are presented in Fig. |l| (bottom panel) and Fig. ^ (mid- 
dle panel). The solid [P(WwS2 + Wd)] and dashed curves 
[P(Wws)l of panel b represent optical-potential calculations 
using the parameters obtained in the ^^O+^^C fit but chang- 
ing the projectile mass. 



FIG. 6. Same caption as for Fig. | but for the ^**0-f ^^C 
system. 

FIG. 7. The position of Airy minima (Oc.m. versus -Ec.m.) 
of the ^^O+^^C system. The solid curves are results of the 
fit to the data using the second-order polynomial as the 
fit function. 

FIG. 8. Same caption as for Fig. § but for the ^*0-H^^C 
system. 



FIG. 9. Real and imaginary parts of the phenomenolog- 
ical optical potentials obtained for the ^^C-|-^^C, ^''O-f^^O, 
^^0-|-^^C, and ^*0-|-^^C systems at three different energies. 
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